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SUMMARY 


The use of a magnetic field to control the motion of electrically 
conducting fluids Is studied. The boundary-layer solutions are found 
for flow over a flat plate -when the magnetic field is fixed relative to 
the plate or to the fluid. The equations are integrated numeric ally for 
the effect of the transverse magnetic field on the velocity and temperature 
profiles, and hence, the skin friction and rate of heat transfer. 

It Is concluded that the skin friction and the heat-transfer rate are 
reduced when the transverse magnetic field Is fixed relative to the plate 
and Increased when fixed relative to the fluid. The total drag Is 
increased in all the cases studied. 


INTRODUCTION 


It has been said that a fluid at a very high temperature Is like a 
universal solvent 'which cannot be contained. A possible method of con- 
taining this fluid Is suggested when It is noted that at such a high 
temperature It would surely contain ions and quite probably also free 
electrons. The fluid would then be an electrical conductor. The invis- 
ible hand of electrical and magnetic fields can then be used to Induce 
forces on the fluid such that It is prevented from coming in direct contact 
with a wail which It would dissolve. 

A somewhat similar technique has been in use for some time In the 
metal purification industry. It employs a high-frequency magnetic field 
which causes eddy currents in a lump of molten metal which In turn react 
with the Imposed magnetic field. The metal is thereby suspended In space 
if the imposed magnetic field Is made strong enough. 

Another example Is the so-called "perhapsatron” described briefly in 
reference 1. A gas In a doughnut shaped container 1 b heated to a high 
temperature by an electrical current discharge. Through the use of the 
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Interaction of the resulting ion current and a magnetic field, the hot 
gas is prevented from coming in contact with the surface of the vessel* 

An application of similar principles was used in thenoonuclear fusion 
experiments in the Soviet Union. The techniques and results described 
very "briefly in reference 2 indicate that it was possible to keep the hot 
fluid from the walls and to concentrate the hot gases quickly so as to 
generate a focusing shock wave. 

These examples indicate the possibilities which may be realized 
through the use of a magnetic -field force on a flowing conducting fluid. 

A number of situations exist in aerodynamics wherein a magnetic or elec- 
trical field might be used to alleviate higgi convective heat-transfer 
rates to a wall. Such problems arise in the flow of air in the boundary 
layer and in the vicinity of stagnation regions of a miBsile or airplane 
moving at very high speeds. If the velocity is high enough the air is 
ionized and, therefore, electrically conducting. See, for example, 
references 3 and * 4 . Other examples are those associated with the flow 
of the combustion products (which are hot and generally electrically 
conducting, refs. 5 through 7) In the propulsion units of aircraft. 

The extent to which a given heating problem can be alleviated is as 
yet difficult to determine. It would be well if a few theoretical results 
were available to evaluate the effect of a magnetic field on the drag and 
heat-transfer rate. An attempt to include ail- the aerodynamic features of 
the flow of air over or inside a missile plus the magneto -hydrodynamic (or 
more correctly, magneto -aerodynamic) effects would render the problem so 
complicated as to make its solution difficult if not impossible. It Is 
felt that a few simple basic theoretical solutions would point out the 
advantages or disadvantages and yield a rough estimate of the various 
quantities entering the problem. It Is the purpose of this report to 
present the results of several such solutions which. It Is hoped, will 
extend the information already available bo that such estimates can be 
made. Several applications of the results are indicated. 
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PRINCIPAL SYMBOLS 

magnetic induction, lines/sq in. 
local skin-friction coefficient 

heat capacity of air at constant pressure, Btu/slug °P 
electric field intensity, volts/in. 

u? 

total energy, CpT + 
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error function. 

complementary error function, 1-erf 

dimensionless Blasius boundary-layer stream function 

dimensionless stream function, equation ( 26 ) 

force, lb 

local heat-transfer coefficient, Btu/sec ft 2 °F 
magnetic intensity, , (ampere ) (turns ) /in . 
electric current, amp 

electric current density in the fluid., amp/sq in. 


electric current density associated, with generation of 
magnetic field., amp/sq In. 

thermal conductivity, Btu/sec ft 2 °F/ft 

characteristic length 

ffB o 2 

magnetic parameter, — — , per in. 
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magnetic Reynolds number, 
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Laplace transform variable 
time, sec 
temperature, °R 

velocity components In x,y,z directions 
Laplace transform of velocity u (see eq. (14)) 

velocity, ft/sec 

coordinate axes, x alined with free stream 
displacement thickness of boundary layer 



excess charge density, coulomb/cu in. 

magnetic viscosity, — ~ 

4rtop 

coefficient of viscosity, slugs/ft sec 

magnetic permeability, (volt) (second) /(ampere) (inch) 

kinematic viscosity, £ , ft 2 /sec 
x 

density of fluid, slugs/cu ft 
conductivity, mhos/in. 
shear Btress, lb/sq in. 
stream function 


Subscripts 


basic quantity 
free stream 
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B 

magnetic induction 

e 

electron 

w 

at plate, tj = y = 0 


Superscripts 

f*) 

vector quantity 


( )', ( )", ( ) ,n derivatives 

HISTORICAL REVIEW OF MAGNETO -HYDRODYNAMICS 


Hie earliest known published, works treating a problem in the flew of 
an electrically conducting fluid through a magnetic field are those of 
TTftWanwrm and of Hartmann and Lazarus in references 8 and 9- Since that 
time a number of theoretical and experimental studies have been carried 
out. Some of these pertain to the flow of a conducting fluid but a larger 
number relate to the dynamics of Ionized clouds or stars. Ho attempt 
will be made here to Include a complete discussion of the latter group. 
However, much can be learned by studying the papers an the flow of con- 
ducting fluids. Haese papers will be divided into several groups according 
to the type of problem treated. 

1. Flow in a channel (refs. 8 through 19)* 

2. Flow about bodies of various shapeB (refs. 20 through 22). 

3. Related papers on astronomy (refs. 23 through 26). 

The various papers will be reviewed briefly on a grot?) basis and reference 
to particular papers given only in special Instances. 


Flow In a Channel 

Following his Invention of the electromagnetic pump 1 in 1918, 
Hartmann, reference 8, developed equations which describe the steady- 
state flow Inside a channel In the presence of a magnetic field or a 

■Electromagnetic pumps are used in atomic energy plants to transfer 
radioactive solutions from one place to another (e.g. , see p. 3 2 of Atomic 
Power - a Scientific American book, Simon and Schuster Inc., New York, 

1955)* 
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magjnetic and. electrical field. Interest was stimulated In problems of 
this type when It was found in reference 9 that a turbulent stream could 
be stabilized to the extent that it could be forced to return to a laminar 
flow. Such a sequence of events is opposite to the usual form of transi- 
tion when a magnetic field is not present. Further investigations by other 
authors, references 10 through 19 , substantiated these results and led to 

the introduction of terms such as "magnetic viscosity," A = — ~ , and 

_ r - 

"magnetic Reynolds number," M = B2. J— or Q = — — , 

Thg stabilizing force exerted on a conducting fluid by a magnetic 
field perpendicular to the flow direction comes about because of the elec- 
trical current which is induced. The portion of the fluid moving most 
rapidly crosses more magnetic lines of force than its neighbor which is 
moving at a slower rate. Since the electrical potential induced locally 
is directly proportional to the number of lines crossed per unit time, a 
potential or voltage difference is brought about between the two lumps of 
fluid. If the fluid is conducting, an electrical current will flow between 
the two points . By virtue of the resulting induced magnetic forces the 
more rapidly moving fluid is decelerated and the Blower moving fluid is 
accelerated. 

Stabilization is also achieved because of the fact that the conducting 
fluid experiences a force which resists the motion across the magnetic 
lines of force. For example, if the main fluid flow direction is parallel 
to the magnetic lines of force, no interaction takes place between fluid 
and magnetic field unless same deviation or instability arises . In this 
event, a restraining force will be induced which is proportional and 
opposite to the velocity component perpendicular to the lines of force, 
tending to damp out the Initial instability. 

Ole types of stability described above were studied in references 13 
and 14. The transition boundaries are given in graphical form for a range 
of the magnetic parameter. It was observed that a considerable Increase 
in laminar run could be achieved at the higher values of the m agnetic 
parameter. 


Flow About Bodies of Various Shapes 


The papers describing the external flow of a conducting fluid about 
bodies are amaH in number. The three references listed, 20 through 22, 
describe the flow of a conducting fluid around a circular cylinder and a 
sphere in the presence of a magnetic field; the third paper treats the 
stab ili ty of the flow between rotating cylinders. Streamlines are shown 
in the first two cases and the stability boundary for the latter case. 
Problems in this category are naturally more difficult than the 
unidirectional channel flow problems discussed in the previous section. 
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Belated Papers on Astronomy 


A need for understanding the motion of electrically conducting fluids 
(ionized clouds) has "been felt by astronomers for some time. Die dimen- 
sions and quantities considered are generally not reproducible in the lab- 
oratory or in the earth’s atmosphere; however, some of the equations and 
results may be applied to aerodynamic problems. 

In reference 23, Batchelor investigated the possibility of spontaneous 
magnetic fields arising in a conducting medium as a result of turbulence. 

He found that unless 4npov > 1 such fields cannot occur. In all foresee- 
able problems in aerodynamics 4upcrv is many orders of magnitude less 
than 1, ruling out spontaneous magnetic fields in a fluid due to 
turbulence. 

Hie manne r and velocity at ■which pressure waves are propagated in 
an ionized fluid and other cosmic problems are treated in references 2k 
through 26. Hie results of these studies an wave propagation may be 
useful in studying interference problems on high-velocity aircraft. 


MAGNETO -AERODYNAMIC EQUATIONS 


Hie equations which describe the flow of an electrically conducting 
fluid in the presence of magnetic and electrical fields will now be dis- 
cussed. Hie equations are simplified to the extent that the flow is 
assumed to be incompressible and to have constant properties (Cp, p, k) 
throughout the flow field. In all cases except one, the conductivity of 
the fluid is assumed to be constant. 


Momentum-Transport Equations 


The momentum-transport equations for the flow of a viscous 
incompressible fluid consist of a combination of the force terms arising 
from the excess charge density 0, and induced magnetic effect due to -Hie 
motion of the conducting fluid through magnetic lines of force (see, e.g., 
refs. 23 end 27) and the usual Navier-Stokes relations (see, e.g., 
ref. 28). In equation form, these force terms are. 


3 ? = 


0E + j X B 


where the term 0E results from the electrostatic force on the excess 
cha rges due to the presence of an imposed electric field if! Hie second 
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term describes the force on the fluid, due to the interaction of the 

electric current, ^ In the fluid, and. the magnetic induction if". The 
differential equation can he written in vector notation as 


EU 

Dt 


— - ^ T X E 
P P 


i grad 


vV 2 U 


(1) 


where p is the fluid density, 9 is the excess charge density, and j I 
is the magnetic permeability. The substantial derivative D/Dt is 
associated with the use of spatial coordinates. A simplified version 
of Maxwell's equations appropriate for the present problem is 


div H =* 0 ; div J = 0 ; pH = B (2a) 

curl E = (2b) 


A reduced form of the generalized Ohm's law 1 b, from reference 29, 


J 


c + PU 


X H) 


(2c) 


The more complete form gLven as equation 2-22 in reference 29 includes 
additional terms accounting for gravity and pressure gradient effects. 
Both of these factors will be assumed negligible in this analysis. It 
will also be assumed that the conductivity, a, does not have directional 
properties. 


The continuity equation is 


div U* = 0 


( 3 ) 


Equation (2b) expresses the relation between the magnetic intensity 

H and the total electric current density J. The distinction between 

the current densities and ? lieB in their location and action. Con- 
sider first the situation when the fluid is stationary and the electric 

field intensity iT is zero. The current density in the fluid is then 

zero. The magnetic field strength Hq or magnetic ind u ction B Q = pH 0 

is a result of an electric current ^ outside of the fluid. Hie 
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current density J Q could be tbou^it of as an electric current in a coil 
or electromagnet outside the flow region producing a magnetic field of a 

given number of lines of force per unit area. The relation between J n 
■> ° 
and EL, is given by equation (2b). The basic magnetic field strength will 

be designated as H Q . 

^ When a conducting fluid moves through the magnetic lines of force, 

B 0 , the positive and negative charges are each accelerated in such* a way 

that their average motion gives rise to an electric current J = crtJ x B, 

where B = B n + b. The quantity b is the magnetic induction resulting 

■» 

from the electric current j in the fluid. In this analysis b will be 
considered as a perturbation on the basic field strength and negligible in 

comparison with B^; that is, from equation (2b), 


curl(B 0 + b) = hitp.(J 0 + j) = itrtyJT 


or 


«> 


curl K, 




■> 

curl b = 


•> 

where B 0 » b. 

If an electric potential E were to be imposed across the flow field, 
the current density would be changed by =* oE*. However, the boundary 

conditions will be taken bo that the impressed potential, E, is zero. 

Therefore, = oif” x that 1 b, the electric current is proportional to 
the voltage generated by the relative motion of the fluid and magnetic 
field. 

Thf> fluid is assumed to be ionized and thereby an electrical 
conductor. However, within any small but finite volume the number of 
particles with positive and negative charges are nearly equal. The total 
excess charge, 9, in this small but finite volume will be taken as equal 
to zero. It has been shown in reference 17 that a very small increase 
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In viscosity is realized due to the Increased ordered-mLgration velocities 
of positive and negative Ions in the presence of the electric field. This 

effect -will he neglected and the electrostatic force term, 0E will he 
assumed to he identically zero. 

The induced electromagnetic force term, J x B 0 will now he related 
to the local velocity vector. The procedure used in references 10, 19, 

and 22 which results in the magnetic parameter (ratio of electro- 

pU 2 


magnetic and inertia forces) and the magnetic viscosity X 


WafT 


will 


not he used. Instead a procedure will he foil .owed which will lead, to 


oB 2 1 
pU 


, used hy references 7, 9# H* 12, 


the magnetic parameter, Q 
and 13. 

A relation between j x B and the velocity Is obtained from 
equation (2). It will he assumed that: 

(l) Magnetic field lines are perpendicular to free -stream 
velocity. 


(2) Permeability {T is constant throughout the fluid. 

(3) The assumptions already made are that the excess 
charge density 0 and imposed electric field 

Intensity ^ are assumed to he zero. 

(4) Induced magnetic field h Is negligible, linearizing 
the induced magnetic force term. 


The first assumption is not highly restrictive since the perpendicular 
component of the magnetic field 1 b the only part contributing a force 
which will change the laminar velocity profile. The other assumptions 
are reasonable physical approximations. 

If these assumptions are used together with equation (2c), the 
Induced magnetic force term becomes 


4 




* 


■» 

i x B 0 


»-2(?xr 0 ) x f 0 


(5) 


* 
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Prom vector relations , 


fxt n 


- § xt 0 ) xf 0 = f [?(%> • to) - tjt a . t)] 


However, It lias teen assumed that B 0 is perpendicular to the free -stream 

velocity Ug,, (assumption (l) ) . In the problems treated, the local velocity 
direction differs a negligible amount from the free-stream. direction so 

that • Tj” = 0 and the linearized magneto -hydrodynamic force term becomes 


4 » 


i ££.*.«$ 


The parameter 2 m is related to the viscous and magnetic Reynolds numbers 
of references 8 through 12 by 


oB 0 2 l 


oBo a l 2 


Equation (l) may then be written as 


+ °^ s> ■■IT + — grad p = vV 2 ^ 

Dt P P 


In the development of equation (7)* it was assumed that the velocity 
of the magnetic field is zero. Since only the Inductive force is being 


2 Prom electric motor equations, it can be shown that 

_ 2 force 

0B o =, — 

l x velocity 


q -^o ^ _ magnetic pressure 
pTJgj, dynamic pressure 


and therefore 
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considered, the second term becomes 



for a magnetic field 


In uniform translation, -where Ug is the -velocity of the magnetic field. 


When electrical currents flow In a plasma it is well to define the 
conditions at the solid boundaries of the flow field. The problems 
treated are assumed to be two-dimensional with the electric current flow- 
ing from and to Infinity parallel to a flat plate and perpendicular to the 
stream direction. The configuration may also be thought to consist of flow 
in a wide channel with the two side wallB as conductors connected by a con- 
ductor of negligible resistance located outside of the Influence of the 
flow field. The flow in the middle of such a channel should approximate 
the two-dimensional type of flow being analyzed. Such phenomena as charge 
accumulation and boundary effects at the walls will be assumed negligible. 


Thermal Energy Transport Equation 


The differential equation describing the relationship between the 
convection and conduction of thermal energy and the work done on an 
electrically conducting fluid in the presence of a magnetic field is 
found by considering the energy entering and leaving an elementary cubic 
box. The result is given by 



The derivation of equation (8), with the exception of the second term on 
the right of the equal sign, can be found, for example, in reference 28. 
The. term on the left of the equal sign takes account of the en thalp y 
transported by the motion of the fluid or by convection. The first term 
on the right represents heat conducted from one element of fluid to the 
next by molecular motion. The expression, J 2 /cr, is the heat added by the 
electrical current produced by the motion of the fluid through magnetic 
lines of force or by an imposed electrical field (see eq. (2c)). The 
re m a inin g terms arise when work 1 b done on the fluid by either pressure 
or shear forces. The symbol D/Dt is the substantial derivative of the 
particular quantity and V 2 denotes the Laplaeian. 
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If -the general expression for the electric current density, as given 
"by equation (2c), is inserted into the expression for the heat added by 
electrical means, the result is 



( 9 ) 


It has been assumed that the electric field E is zero and that the 
magnetic field lines are perpendicular to the free-stream direction. 

To a good approximation in the problems to be considered, equation (9) 
can then be written as 


^ - OU^O 2 (10) 

where U 2 must be written as (U - IFg) 2 when relative motion exists 
between the magnetic field and the plate. Equation (8) is then 



IMPULSIVE MDTIOff OF A FLAT ELATE 


The impulsive motion of an infinite flat plate 3 in a viscous fluid 
serves as a model for the boundary layer an a semi -infinite fla t plate. 

The advantages of working the problem of the impulsive motion of a flat 
plate are that it 1 b simple enough to yield a result in closed form and 
suggests the proper choice of parameters to be used on the more complicated 
problem, of a send -infinite flat plate. 

The velocity profiles will be found for the two cases. In the first 
it will be as Burned that the magnetic lines of force are f ixe d relative to 
the plate, and the second relative to the fluid. In both cases the fluid 
is of uniform density and viscosity, and has the same electrical conduc- 
tivity throug h out. A third case, in which the magnetic field is fixed on 

g 

This problem is often referred to as the "Rayleigh, problem" when the 
fluid is incompressible and nonconducting. 
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the plate and a compensating pressure gradient exists in the fluid, will 
he shown to he equivalent to the case wherein the magnetic field is fixed 
relative to the fluid. Hie problem wherein the fluid conductivity varies 
will not he considered. The fluid and plate will he assumed to he ini- 
tially at rest. At time t = 0, the plate will move impulsively with a 
given velocity. 


Magnetic Field Fixed Relative to the Plate 



At time t < 0, the fluid, plate, and 
magnetic field are assumed to he everywhere 
stationary. At time t = 0 and for all later 
times the plate and magnetic field are moving 
at velocity u = u^ (sketch (a)). The prob- 
lem is to find the velocity-time history of 
the fluid. 

The velocity in a given y = constant 
plane does not change with x. It can then 
he reaaoned that the vertical, v, and trans- 
verse, w, velocities are zero or negligible 
together with the pressure gradients in all 
directions. Equation ( 7 ) then reduces to 


St P 



( 12 ) 


Since the magnetic field is moving and the fluid is Initially at rest the 
relative motion must he accounted for. Equation (12) is then 


du oB 0 g 

at + P 


(*-«»> 



( 13 ) 


The boundary conditions sire 

u = , y = 0 , t > 0 

u = 0 , y > 0 , t =0 


The Laplace transform of the velocity 


u 


is defined as 
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XL = 



e" st u dt 


(l4) 


Applying the Laplace transformation (see, 
the first term In equation (13) gives 


e.g. , refs. 30 and 31) to 



/ -st chi -j. -st 

e dt = +ue 


00 

+ 

o 



dt 


a bu 


If the other terms of equation (13) are treated similarly, the transformed 
equation is 


dfu 

su = v — — + - 

dy 2 8 




or 


d^tT _ _ / s + V v _ *430 
dy 2 \ v / vb 


( 15 ) 


Hie solution to equation (15) la the sum of the solution of the homogeneous 
equation plus the particular solution, that is. 


* 1^00 




s(s + n^) s(s + n^) 


O x (s)e ** v 



The c onstan t C 2 is chosen equal to zero to fit the boundary condition 
that u must be finite at y = 00. The Integration constant C 1 (s) is 
found from the boundary condition at y = 0; that is, on the plate 


^=0 = ^00 



The Laplace transform of u^ is u^/s. Therefore 


C3.Cs) 


s 

% 


and 

U*, ~ y J ~V 2: ‘ 

u => — -i-2 — + — — — e ^ 
s(s +103.) s + m^ 


(16) 



Sketch (b) 


Equation (l6) can he Inverted 
with the aid of reference 30. 
Combining several terms gives 



( 17 ) 


The symbol erf denotes the 
error function of the argument. 
The error function is discussed 
in reference 31 and tabulated 
extensively in reference 32. 
Typical velocity profiles are 
shown in sketch (h). If m3, 
is set equal to zero, equa- 
tion (17) reduces to the result 
for the Rayleigh problem; that 
is. 





erf 



The velocity is a function of y/ 2 «/vt only and therefore a single profile 
suffices. However, when a magnetic field is acting on the fluid, the 
velocity profiles are not similar because they change with time according 

to The fluid at infinity is accelerated by the magnetic field so 

that the entire masB of fluid is accelerated by the magnetic field. It 
is, in fact, accelerated more rapidly than yhen the only force is the 
viscous action between layers. 
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The skin -friction coefficient is expressed as 



so that, from equation (17) , 

o f (1 

Values of c f are shown in graphical form in sketch (c). When is 

zero, equation (19a) reduces to the skin friction for the nonelectro - 
magnetic Rayleigh problem. 



J 1 I L 



18 


HACA TN 397a. 


The change in the velocity profile is a result of the force exerted 
"by the magnetic field on the fluid. Hie reaction on the unit gene r ating 
the magnetic field is expressed as 


F/unit area = / trB o 2 (u D0 - u)dy = oB 0 2 e“ mit / erf — Z — fly (ijsfo) 


The force per unit area goes to infinity as the upper limit of integration 
goes to infini ty "because the magnetic field extends undimini shed an 
infinite distance from the plate. Thereby, an infinite amount of fluid 1 b 
accelerated, resulting in an infinite force. 


Magnetic Field Fixed Relative to the Fluid 


At all times less than zero the fluid, magnetic field, and plate are 
assumed to he at rest. At time t = 0 the plate begins moving with veloc- 
ity u = Uoo but the magnetic field remains at rest. The differential 
equation is the same as equation (12) because there 1 b no relative motion 
between the fluid at y => °° and the magnetic field. 


— + m.u = 
St ^ L 


Sju 

Sy 2 


( 20 ) 


The boundary conditions ere 



at y = 0 

t > 0 

u = 0 

at y > 0 

t - 0 


Applying the Laplace transformation to equation ( 20 ), eus was done in 
the previous section, yields the transformed equation 


V ^2 = u ( s + 


The solution to this ordinary differential equation is 
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air 


* 


tt = 



n +y, 

+ C 2 e 



■where the constant Cg is set equal to zero "because of the requirement 
of a finite velocity at y = oo. The integration constant, C^Cs), is 
found from the "boundary condition at y =» 0, t > 0, as 

- T 

Equation (20) can "be written as 


u = 


(s + m^) -n^ 



This equation can be inverted with the aid of reference 30 as 



The symbol erfc denotes the 
complementary error function 
defined as 1 - erf (see, 
e.g. , P* 370 of ref* 31). 

The velocity is once again 
dependent on more than one 
parameter so that a single 
similar profile cannot be 
drawn. Several profiles are 
shown in sketch (d). The 
fluid at y = 00 is not dis- 
turbed in this case because 
the velocity acrosB the 
magnetic lines is zero. In 
the vicinity of the plate 
the induced magnetic force 
counteracts the accelera- 
tion force of viscosity 
resulting in an increased 
rate of shear at the walls 
as expressed by the skin- 
friction coefficient 



0 a jb is 



Sketch (d) 
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Even at time, t *= », a friction force acts on the plate - a situation 
contrary to the previous case and the Rayleigh problem. Values computed 
by equation (22) are shown in sketch (c). 

The force on the magnetic field is given by 


F/unlt area 



cjB 0 2 u dy 


(23) 


where u is given by equation (21). The force is finite in this case. 

Another case which may be of interest arises when the magnetic field 
is fixed on the plate but a pressure gradient compensates for the action 
of the magnetic field at y - 00 . A flow of this character may be imagined 
to exist in the early stages of boundary-layer formation on the walls of 
a two-dimensional channel. The conducting fluid 1 b assumed to be pumped 
through a channel containing a stationary transverse magnetic field. A 
pressure gradient arises automatically to compensate for the resistance 
of the magnetic field. The equations then describe the boundary layer as 
a function of time, after the pump has started. For this case equation (7) 
reduces to 


It + 


O + 


i & n. v 

P dx 


The pressure gradient is 

mentioned above at y = 00 . At 


found by examining the boundary condition 
y = 00 


du = ^u _ 
dt dy 2 


Therefore, at y « <», 


+ T 


& 

5x 


0 


and the above equation becomes equation (12 ). Since the boundary 
conditions are the same as when the magnetic -field was fixed relative to 
the fluid, the results shown in sketches (c) and (d) apply here also. 
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BOUNDARY LAYER ON A SEMI-INFIHIXE FLAX PLATE 


la the analysis of the boundary layer on a flat plate moving 
impulsively, it was found that the shape of the boundary-layer velocity 
prof ile changes with time, necessitating a new calculation at each instant. 
Therefore, a simple transformation of the coordinates would not yield a 
similar solution wherein all profiles could be drawn as one. It is then 
to be expected that in the flow over a send. -infinite flat plate in the 
presence of a magnetic field, a si mil a-r solution will not exist. Since 
the boundary-layer solution of Blasius is not known in closed form and 
the velocity profiles at two different stations along the plate are not 
similar, a numerical solution to a series expansion is found by the method 
used in reference 33 . The first- or zero-order term is found to be the 
same as the Blasius solution which is well known and is tabulated, for 
example, on page 103 of reference 28. The first nonzero term will be 
calculated for the cases studied here. 

The boundary-layer velocity and temperature profiles will be found 
for three cases. The conductivity of the fluid is assumed to be uniform 
throughout the flow field in the first and second cases. A particular 
variation is assumed in the third case. The magnetic field is fixed 
r el.at l.TC to the plate in the first and third cases and relative to the 
fluid in the second case. A fourth case is shown to be equivalent to 
the second. 


Magnetic Field Fixed Relative to the Plate 


Velocity profile. - The magnetic field 
lines of force are assumed to be perpendicular 
to the free -stream direction and to begin at 
the leading edge of the plate as Indicated in 
sketch (e). An experimental setup may be 
imagined to consist of a wing or flat plate 
in a stream of conducting fluid such as 
mercury, sodium, or salt water. A magnetic 
field is impressed aero Be the entire channel 
using either a permanent or electromagnet. 



Sketch (e) 


On the basis of the usual boundary-layer assumptions, equation ( 7 ) 
reduces to 


du Su 

u — + v — + m x u 
ox oy 


a 2 u 

dy 2 


(24) 


The boundary conditions are that u = 0 at y = 0 and u = u^ upstream 
of the plate. At y = °°, du/dy = v = 0 and chi/dx = -n^. Equation (24) 
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is the s ame as the incompressible boundary-layer equation except for the 
addi ti onal term n^u. As was done in the Blasius boundary-layer solution, 
the transformation, 


I 

il 



(25) 


is in troduced. Hie stream function, i(r, is defined as 





3 /2 

+ ^mx f + mxf 2 + (mx) f £ 


+ (mx) E f 


4 


■] 


( 26 ) 


where the functions f Q , f x , f 2 , . . . are functions of i\ only. The 
series expansion in ,/mx is suggested by equation (21), where the exact 
form of m is yet to be determined. 

The velocity components and their derivatives are found from 
equations ( 25 ) and ( 26 ) as 


dy d£ dy dr) dy 


u = u^f 0 * +^x f x ' + mxf 2 ' + (mxf^fg 1 + (mxff/ . . .] 


(27) 


v = _ = _ <*fc 

dx dg dx dt] dx 


T-a 

2 V X 


f 0 ' + ,yix f x ’ + mxf 2 ' + (mx) 3 ^* + (mx) a f 4 ’ + 




- J ^ ^f 0 +2^/^ f x + 3mxf 2 + 4(mx) 3/2 f 3 + 5(m*) 2 f 4 + • • •] ( 28 ) 
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Similarly, 

_ JH 2s [: f c u + fl n + mxf 2 M + (mx) 3/2 f 3 n + faff/ + ...] + 
ax 2 X |_ J 

^ £/mx f 1 * + 2axf z ' + 3 (jnx) 3/a f 3 ’ + Mmx^f/ + . . .J ( 29 a) 


a&.T^ / 3 S 

dy ^ vx 


f 0 w +^x f x " + mxf a " + (mx) 3/a f 3 n + (mx) 2 ^" + . . .J 


(29b) 


[f 0 «* + f x *» + mxf £ «* + Gax) 3/2 f 3 ‘» + faff 4 »* + . . .] 


( 29 c) 


When equations ( 27 ) through ( 29 ) are substituted Into equation (24) and 
the like powers of 


mx = 


cBo 2 _ if 


P^» Re x 

are equated, the following differential equations result: 


Of «• 

Ci o 


2f x ,n 


2f a m 


2f 3 ,n 


2f 4 m 


(30) 


» -f 0 "f 0 (31a) 

- VV " f o f x" - 2f X f o" (31b) 

- 2f 0 ’f a * + f x *f x ‘ - f 0 f £ " - f x f x H - 3f 2 fo" + 2f 0 ‘ ( 31 c) 

= 3W + 3f 2 ,f x* " f o*s" - 2f x f a n ~ 3f a f x" “ kf a*o" + 2f x (31*) 

= 4f 0 *f^ + 3W + 2f 2 ' f a 1 + £ 3 'f x ' - £ 0 £/ - 2f x f 3 ” - 
3£ 2 f 2 " - 4f 3 f x " - 5f 4 f Q n + 2fg* (3le) 
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The boundary conditions are 

f 0 = f x = f a = f 3 = f 4 = . . . = 0 at t[ = 0 

f o * = f x * = f 2 * - f 3 * = f 4 * = . . . = 0 at Tj = 0 

f 0 * = 1 and f 2 * = -1 at T) = « 

f X ‘ - f 3 ‘ = f 4 * = fg' = . . . a 0 at TJ = CD 


From the boundary conditions and equations (31b) and ( 31d) , the solution 
to fx and f 3 is taken to be zero throughout the flow field. This fol- 


lows from the fact that fx and f 3 , . . . are zero on all the boundaries 


and the differential equations are satisfied if f x = f 3 = 0 . 
if the solution is unique, it must be zero (eigen solutions are 


Therefore, 
excluded) . 


Likewise, it is seen that the solutions for f 2 , f 4 , . . . are not zero. 
The differential equations may be written as 


f 2 ,n = fo'fs' - | f 0 f 2" “ | faf 0 " + fo’ (32a) 

U w = 2f 0 , f4* + f 2 'f 2 * - | fof4 B - | fsfa" - | f4fo" + fa’ (32 b) 


The differential equations for f 2 and f* are linear ordinary 
differential equations with variable coefficients. Since the coefficients 
(Blaslus solution) are not known in closed form, the equations must be 
Integrated numerically. Only the function f 2 of the series will be 
found because the work involved is sizable and it 1 b the intent of this 
study to find the over-all trends and gross effects rather than precise 
results . 

The Runge-Kutta method was used to integrate both the homogeneous 
and nonhomo geneous forms of equation (32a). The integration is started 
at the surface of the plate using two boundary conditions and then assum- 
ing a value for £ z ". The correct solution is found by a proper combina- 
tion of the numerical solutions to the homogeneous and nonhomo geneous 
equations. The factor by which the homogeneous equation is to be multi- 
plied is found from the boundary condition at q = ». The numerical 
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results for fa are tabulated 
In table l(a) and velocity pro- 
files for several values of mx 
are shown In sketch (f). Hie 
velocity profile for mx = 0.5 
Is drawn to Indicate trends and 
is probably not precise since 
higher order terms in equa- 
tion ( 27 ) may have became 
notable. 

Hie skin -friction coeffi- 
cient and the displacement 
thickness are 


c f - 


= 0.664 - 1.788 mx + . . . 

# ie x 


(33a) 


8* = mxy+ (l. 73 + 0. 54 mx) /— 

^00 

(33*) 



The drag on the plate which arises because of the magnetic -field- 
fluid Interaction is 


F/unit area 



y->co 


aB Q a u dy 


Substituting equation ( 27 ), 




F/unit area = oB 0 2 ^ u^vx / (f 0 * + mxf 2 ’ + . . . )dTj 

o 

= dB 0 Z (f Q + mxf 2 + . . . ) o n_>> “ (33c) 

Hie force is finite only for a conducting layer of finite thickness. 

Temperature profile .- When the pressure gradients are reasoned to be 
everywhere zero and the usual two-dimensional boundary-layer assumptions 
are made, equation (ll) reduces to 
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St 


St 


puCp g + pvCp ^ 


Sy 


OBrTu' 


2„2 


£e^£t . 

Pr Sy 2 



(34) 


If total energy is defined as 


E 



(35) 


equations (24) and (34) can "be combined to yield a differential equation 

for the transport of total energy Ej that is, if equation (24) is multi- 
plied by u, added to equation (34), and the Prandtl number, Pr, is assumed 
to be 1.0, 


u 


SE SE 

Sx Sy 



(36) 


Although the temperature distribution is altered by the presence of the 
magnetic field because u and v are affected, the total energy of the 
conducting fluid is not. The kinetic energy removed by the force of the 
magnetic field is exactly equal to the heat generated by the electric 
current, independent of the field strength Bq and the conductivity cr. 

An analogy may be drawn by considering a similar situation for an elec- 
tric motor. Assume the motor to be rotating at a given speed. The power 
is then turned off and the armature is short-circuited. The motor will 
decelerate at a given rate depending on the field strength, resistance of 
the circuit, and rotational speed, Just as in the fluid flow problem. The 
temperature rise of the wires, insulation, etc., however, depends only on 
the initial kinetic energy of the rotor. If the conductivity of the wires 
(or fluid) is reduced, the time required to stop the motion will be 
increased but the total j 2 /cr or joule heat is not changed. 

A first-order es tima te of the influence of the magnetic field on the 
temperature profile can be found from equation (34) together with equa- 
tions (25) through ( 29 ) and table l(a). Assume that 


T = T q (ti) + T 2 (q)mx + T 4 (Ti)nft^ + . . . (37) 


where, T a , T 4 , 


are functions only of tj = y 


H22 

vx 


and 


T 


To = T. 


. . = 0 as in the part of the problem dealing with the 
velocity profile. The derivatives of the temperature with respect to 
x and y are found as 
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ST 

St 

M + 

St St^ 

Sx * 

' Sg 

Sx 

St] Sx 

St s 

3 1 

mxT 

+ 2m^x^T 

Sx 

X 

2 . 


St _ 

St 

^21 = 


dy 

Stj 

Sy 

^vx v ° 


. - — (T 0 * + wsZ z x + m^x 2 ^* + 



When the various expressions are inserted into equation (3*0 &nd teims 
containing nice powers of mx are equated, the differential equations 
for T 0 and Tg are 


Pr 


+ 




(39) 


V + Tg£o 

Pr 2 


- fo% + 1 w + 



Of “f 

^ r O X z 


= 0 


(to) 


The boundary conditions are 



= Sir 

at 

H 

- 0 

To 

= 

at 

11 

=a 00 

m 11 
J -2 

= Tg l = 0 

at 

11 

= 00 

T 2 

= 0 

at 

n 

= 0 


The solution to equation ( 39 ) was first found by Pohlhausen and is 
discussed on page 246 ff . of reference 28. The function T 0 (tj) is the 
temperature profile for a boundary layer on a flat plate and is given 
by (assuming Pr = 1.0) 


T 

x o 


= T + 
-*•00 T 





( 41 ) 
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The derivative of T Q to "be inserted into equation (40) becomes 


T 


t 

o 



(42) 


For simplicity, the temperature of the plate _ will be assumed to be the 
same as that of the fluid far from the plate; that is, T w = T M . 
Equation ( 40 ) then becomes (Pr = 1.0) 



f o " 


fo'*a 



f f 11 
X 2 A 0 


2Cp 




Equation (43) wbb integrated 
numerically by the same method 
used for equation ( 32 a) and 
the results are tabulated in 
table 1(b). Several tempera- 
ture profiles are shown in 
sketch (g). The quantity of 
heat transferred to the plate 
per unit time is 



(44a) 


The local convective heat- 
transfer rate is 


Sketch (g) 
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h 5 M*£r. (o. 664 - 0.704 rax -...+.. .) (44b) 

u^/aCp 


la the limit as the conductivity or field strength vanishes (m — ^ 0), 
equation (44b) reduces to 



k£oXP 


(45a) 


or by virtue of equation (33a) 



1 Cp 

i c A^ 


(45b) 


Equation (45b) is often identified as Reynolds analogy. A similar relation 
of equal simplicity was not found for the problems treated in this paper. 

Hie differential equation for flow over an insulated plate is obtained 
from equation (40) by introducing the proper expression for T 0 ' . Hie 
solution of equation (39) for the insulated plate case is (see, e.g. , 
ref. 28) 



(46a) 


(46b) 


Equation (40) then becomes, for the insulated plate case (Pr = 1.0), 


qi If Ot, 2 

rp n , 2 O j. »m O 13 *” -P p 1 -P 

i 2 + 2 r O 1 2 r 2 r ° Zr 


« + 13,00 f -p * 2 + 2f n f - 
+ Cp Vo + 2r o r 2 )- 


(46c) 
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Die temperature profiles 
(table 1(c)) for several values 
of mx are shown in sketch (h). 
It is to be noted that the 
recovery temperature is not 
altered by the magnetic field. 
Such a result is to be expected 
since the magnetic field and 
the fluid at the wall do not 
interact. 


Magnetic Field Fixed Relative 
to the Fluid 


Velocity profile .- In this 
case the magnetic lines of force 
move past the plate at the free- 
stresm velocity. Die fluid, 
■which is decelerated by the 
viscous action of the plate, 
receives a push from the mag- 
netic field which counteracts 
the viscous effect. Such a 
situation exists when a plate, 
vane, or wing moves through a 
stationary electrically conducting fluid and magnetic field. Die bound- 
ary layer on a wing of an aircraft flying in the ionosphere over the mag- 
netic poles of the earth could also be considered to correspond to this 
case. Under these circumstances, equation ( 7 ) reduces to 



Sketch (h) 


(u-O-vgji 


If the transformation of coordinates and the stream function described by 
equations ( 25 ) and ( 2 6 ) are introduced into equation (47), a set of ordi- 
nary linear differential equations with variable coefficients is found for 
f x , f 2 , f 3 , f 4 , . . . Die boundary conditions are 


f Q a a ^2 3 ^3 3 • • • = 0 at n = 0 

V = = f 2 ' =■ fa* = • . . = 0 at T| = 0 

f o ” 1 and fx* = fj 2 * = f 3 1 =» . . . = 0 at T] = 00 
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As In the previous case, it is found that f x = f 3 = . . . = 0. OIhe 
differential equations for f a and f 4 are t h e n 


f 2 '" = fo’fs* - | fo’fa" “ | fsfo" + V - 1 (48a) 

f 4 «* = 2f 0 'f4' + fz'fs? - \ f 0 f4 U - | fsfs" - | f-tfo" + (48b) 


For the reasons pointed out in 
the previous case, only the 
solution for the function 1 2 
was found. -OIhe Integration was 
carried out by the same method 
described in the previous sec- 
tion and the results are tabu- 
lated in table Il(a) . Several 
velocity profiles are shown in 
sketch (i) . The local skin- 
friction coefficient is 


0.664 + 2.293 mac + . . . 

(49a) 


and the displacement thickness 
Is 



B* = ( 1.73 - 2.21 mx 



(49b) 


The disturbance caused by vIbcous action on the plate results in a 
force on the magnetic field, •which is given by 


ny->°° 

F/unit area = / cffl 0 (u^, - u)dy 
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From equations ( 27 ), (66) , and. table Il(a), 


F/unit area 



- 2 .9 0B o 2 ^ u^vx 


•)4ti I 


(49c) 


The force would appear as a reaction on the unit generating the magnetic 
field and not on the model which disturbs the fluid. 

Equation (47) is identical with the differential equation for the 
case when the magnetic field is fixed with respect to the plate and a 
pressure gradient exists Just strong enough to compensate for the elec- 
tromagnetic effects at tj «= ». The boundary conditions and, hence, also 
the numerical results of table Il(a) and sketch (i) can be applied to 
this problem. 

Temperature profile . - Since the magnetic field lines are now moving 
relative to the plate, the expression for the electric current in the 
fluid is 


J - dB 0 (u - O 


The equation for the transport of thermal energy in the boundary layer 
then becomes 




(50) 


If equations (47) and ( 50 ) are combined to describe the transport of total 
energy in the boundary layer (Pr =» 1 . 0 ), 

pu ~ + pv + u t30 (u - u^JoB 2 = |i ( 51 ) 

dx dy dy 


where, as before, E =» CpT + — . It is seen that the presence of the 

magnetic field adds energy to the local element by doing work on the 
fluid retarded by viscosity. The magnitude of the work depends on the 
velocity, electrical conductivity, and field strength. 
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When the transformation of coordinate s , equation ( 25 ) , and the series 
expansion, equation (37) # are inserted into equation ( 51 )* a set of dif- 
ferential equations is obtained for T 0 , T 2 , T4, . . . The relation for 
T 0 is the same as equation ( 39 ) and that for T 2 is 


T P ” f 0 3 u» 2 r, 

+ trL {1 ’ 


+ 2f "f. 


■I. 


I \ 


When the expression (42) is inserted into equation (52a) and it is assumed 
that T w = ^ and Pr = 1.0, the result is 


Vo"! 


- - f , 


( 52 b) 


The boundary conditions are 



- 0 

at 

T) - 0 


- 0 

at 

8 

1! 

P 


The numerical integration 
of equation ( 52 b) "was carried 
out by the method used earlier 
and the results are tabulated 
in table II (b). The temperature 
profiles shown in sketch (j) 
correspond to the same values 
of mx shown for the velocity 
profiles in sketch (i). The 
temperature gradient and the 
heat transfer to the wall are, 
respectively. 



Sketch (3) 
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2 /— 

£si- - 2° (0.332 + 2.022 mx + 

2Gp Vvx 


• •) 


a E^s££- (0.664 + 4.044 mac + . . .) 

ti 00 2 /2C p 2^Re x 


(53) 


The recovery temperature for an Insulated plate Is altered by the 
relative motion between the magnetic field and the fluid at the plate. 

Th e differential equation for the perturbation profiles are obtained by 
inserting equation (46b) into equation ( 52 a)* 

T z " + V - £ 0 'T Z - f a f «f « + ^ [(1 - f 0 ') 2 + 2f 0 ,, f s ,, l - 0 

2 2Cp °P 4 J 

(54a) 

Several typical temperature profiles (table 11(c)) are shown in sketch (k). 
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The recovery temperature at •Hie wall is changed due to the J 2 /a or Joule 
heating caused by the relative motion between the magnetic field and the 
fluid at the plate; that is, the recovery temperature at the wa,ll is given 
by 

T r = ^o + l*6l4 ^ sa - mx . 

°P 

2 

T r = + — (1 + 3.228 mx .) (3*fb) 

2Cp 


The heat transferred to the plate may also be written as 
q = h(Tr - Tv) > where T r is the recovery temperature given by equa- 
tion (54b) and T w is the wall temperature. Since the magnetic field 
increases T r as well as t, the heat-transfer rate 1 b increased in 
two ways. 


Magnetic Field Fixed Relative to Plate - Fluid of 
Variable Conductivity 


Velocity profile .- numerous problems in fluid dynamics exist wherein 
the electrical conductivity of the fluid or the strength of the magnetic 
field c h anges with distance from tie surface. The magnetic field strength 
can be varied somewhat by the design of the apparatus. It is possible 
to adjust not only the strength of the field but also the variation with 
distance from or along the surface within certain physical limits. A 
first approximation to the shape and magnitude of most magnetic fields 
would be a uniform magnetic field of constant magnitude throughout the 
flow field. Such a magnetic system will once again be chosen for this 
case. 


The electrical conductivity of the fluid would be expected to vary 
in a complicated manner in directions along and perpendicular to the sur- 
face and from one situation to the next. The constant value chosen for 
the conductivity in the two previous cases is not realistic. The con- 
ductivity will vary under actual conditions of flow in a boundary layer. 
This will surely change the skin friction and heat transfer from that 
found in the last two cases. A logical variation will be studied in this 
section to see how large an effect a variable conductivity will have an 
C£ and h. 
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It is to be expected that the conductivity will vanish at the outer 
of the boundary layer and reach a ma-iriTniTni somewhere near the surface. 

The equations which can be used 
to determine the number of ions 
(and hence conductivity) in a 
stream at high temperature or 
under circumstances wherein 
combustion is taking place, are 
cumbersome. Therefore, a sim- 
ple, approximate, analytical 
relation for the electrical con- 
ductivity as a function of tem- 
perature between the plate and 
outer edge of the boundary layer 
was not found. The conductivity 
of air after the passage of a 
normal shock wave "has been 
measured in a shock tube and 
reported on page 339 of refer- 
ence 3. The conductivity is 
reproduced in sketch (l). At 
Mach numbers above 20, the data 
are approximated by a straight 
line. 



In view of the lack of knowledge of the variation of a in the 
boundary layer, a linear relation between the velocity decrement and the 
conductivity will be chosen for this case. Therefore, it will be assumed 
that 


a = 



(55) 


where cr 0 is the conductivity at the surface, where u = 0. Whether 
equation (55) Is realistic is uncertain. Inserting the expression (55) 
into equation (24) gives 


u^ + v^+SLSi^-u) 
dx dy Uoo 


- S £ u 

S3 V — 

dy® 


(56) 


As before, an ordinary differential equation for each of the functions 
f x , f 2 , fa, f 4 , f . . . results when the variables in equations (25) 
and (26) are inserted into equation (56) and like powers of mx are 
equated. The boundary conditions are 
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fr 


f - f 2 = f 3 = f 4 


V - fx* - f 2 * - f 3 * = = • • - = 0 

f Q ' = 1 , f x * = f 2 * = f 3 ’ - f 4 ‘ = . . . = 0 when 


when tj = 0 

when T) => 0 


Once f 1 = f 3 = . . . = 0 throughout the flov field and the 

differential equations for f 2 and f 4 become, 

f 2 m = f 0 *f 2 ‘ - | f 2 f 0 n - | f 0 f 2 " + f 0 '(l “ V) 

f 4 ‘" = 2f 0 *f 4 *+f 2 ‘f a « -i f 0 f4 ,n -§ f 2 f 2 n -| W + W+fa’d 


(57a) 

to') 

(57b) 


The numerical Integration 
of equation (57a) was carried 
out by the method used In the 
other cases and the results 
are tabulated In table IXT(a). 
Typical boundary-layer veloc- 
ity profiles are shown In 
sketch (m). He local skin- 
frlctlon coefficient is 

= 0.66k - 0.627 mac ± . . . 


*y Re x 


( 58 a) 


and the displacement thickness 
Is 


8* = (1.73 + 0.97 “3C + . . .) 

v ^00 
(58b) 



Thf» force on the magnetic field is expressed by 


r^V^eo 

F = 8. 85x10" 16 0B o 2 ^ u^vx / (1- f 0 ‘ -nrs£ z l - ...)(f 0 ’ +“3cf 2 r — )d.Ti/sq. in. 
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The Integral Is evaluated through the use of equation (57a) and 
table m(a) as 


F = 0. 738x10 (0.66k + 0.6k9 mx + . . . )lh/sq In. (59a) 


or, integrating from the leading edge to x. 


F =» 0 . 492xKr 16 xcr 0 B 0 2 J u^vx (0.664 + O.390 mx + . . .)lb/ln. width 

(59b) 

Equation (59b) expresses a drag force on the plate which is to be added 
to the friction drag. In coefficient form, the local drag on the plate 
due to the magnetic field is 


Cg = -pE= (1.328 + 1.299 mx + . . .) (60a) 

*/®®x 

When equations ( 58 a) and (60a) are added, the net local force coefficient 
to a first order in mx is 


Cy, = —= (0.664 + 0.708 mx + . . .) (60b) 


The interference of the magnetic field reduces the skin friction at the 
expense of a higher drag. 


Temperature profile . - The thermal-energy transport equation is 
reduced in the manner used in the previous cases. The electrical or 
Joule heating is now described by " 




u 


Equation (ll) is then 


PV 




gpBo 2 


(u» - uju 2 


Fr Sy 2 





2 


(61) 
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If equation ( 56 ) is multiplied 'by u and added to equation ( 6 l), 
the resulting equation for the transport of total energy is identical 
with equation ( 36 ). Once again, the temperature distribution is changed 
but the total energy of the fluid is not altered by the presence of a 
magnetic field. 

The change of variables g, tj, the stream function, if, and temperature 
expansion given by equations ( 25 ), ( 26 ), and (37), respectively, are now 
introduced into equation (6l). A set of linear ordinary differential 
equations for T 0 , T 2 , T*, ... is obtained for this case when like powers 
of mx are equated. Hie expression for T Q is once again the Barne as 
equation (39). The differential equation for T 2 is 

[(l - f 0 *)V 2 + 2f 0 ,, f 2 n ] = 0 

(62) 

If it is assumed that Pr =■ 1 and the expression for T 0 * is inserted 
into equation ( 62 ), the result is 


T. 


2 

Pr 


f _T 1 

O 2 


-=-+ ~~~~ - fo’Ts 


+ - f T 1 
2 2 0 


T3oo 


V + ^r~ - f °’^ + ^ [| - v) + (1 - V )f c 


2f 0 "f £ 


n 


Numerical integration of equa- 
tion ( 63 ) was carried out by the 
method used for the previous 
canes. The results are presented 
in table Hl(b) and several tem- 
perature profiles are shown in 
sketch (n). The coefficient for 
the convective heat-transfer rate 
is given by 




.+. . . ) 

(64) 


( 63 ) 



Sketch (n) 
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The influence of the magnetic field on the heat-transfer rate and 
flTHn friction is not so large in this case as in either of the two pre- 
vious cases (eqs. (33a), (44b), (49a), and (53) ) - these results, 

it is apparent that the ^’hangy* in the heat-transfer rate and the skin 
friction is strongly dependent on the manner in which the conductivity 
(or magnetic field strength) vary through the boundary layer. A magnetic 
field which changes in strength with distance from the surface and a 
variable conductivity are to be expected in the real physical situation. 

tk.* differential equation for flow over an insulated plate is obtained 
by substitution of equation (46b) into equation ( 62 ). 

la” + ^2^1 . f 0 'T z . f 2 f 0 'f 0 " + ^ [(1 - f 0 *)V Z + 2*0%”] - 0 



Sketch (o) 


Thp. temperature profiles for 
several values of mx are 
shown in Bketch (o) and tabu- 
lated in table IH(c) . As 
would be expected, the recovery 
temperature at the wall is not 
altered by the magnetic field 
because there is no relative 
notion between the fluid at the 
wall and the magnetic field. 


MAGNETIC PARAMETER 


Units 


The importance of the 
magneto -hydrodynamic effect 
depends on the size of the 

OBq 2 

magnetic parameter m . 

The usual engineering units 
are: a in mhos /inch; B 0 in 

lines/square inch; p in 
slugs/cubic foot; and u^ in 


feet/second. If each of these quantities is expressed as such, the param- 
eter mx will be dimensionless only if multiplied by a factor which 
relate b the. various dimens i onal quantities. When the correct conversion 
factors are inserted, reference 27 , 
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m = 1. 529x10 “ — - — per Inch 

'(■aws) 


41 


( 65 a) 


or 


m = 6.432x10 


-10 




per Inch. 


( 65 b) 


■where p Q => 0.002378 slugs /ft 3 is the air density at sea level. 

Another method of finding the conversion factor for m is to combine 
the electric motor equations (ref. 34) 


' 

F = 8.8907XU0‘' b B 0 Zi lb 

( 66 a) 

- 

E = BlvxlQ-* voltB =« —■ 
° al 

( 66 b) 


so that 


For " = 8.85O7XLO" 10 ctB 0 2 B3C 

ul 3 in . 4 

or 

m = >1 = Ca.a gTit ms)^ 1 : to ob 

p 0 (p/p 0 ) 1 l E O.OCE378 (p/p 0 )u 


iq . 

= 6.432x10 per inch 

(p/p 0 Ko 


In words, the parameter m expresses the ratio of electromagnetic force 
density to dynamic pressure. 
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The value of m can eaBlly vary over many orders of magnitude . 
However, the density p and. velocity u associated with a particular 
vehicle are generally fixed within relatively narrow limits. It remains, 
then, to adjust the conductivity and magnetic field strength so that tfye 
desired characteristics are achieved. The limits for the magnetic field 
strength lie between zero and a maximum (die counting astrophysical possi- 
bilities of about 10 5 lines/sq in. The electrical conductivity depends 
on the fluid considered and can vary over wide limits from very small 
values for insulators to approximately 10 s mhos/inch for metals. The 
results of several investigations of the conductivity of gases will be 
discussed in the next section. 


Electrical Conductivity 


The information on the electrical conductivity 4 of gases under the 
special conditions of this report is not large. However, it 1 b known 
that the air around a missile will be a conductor if: (l) the energy 
imparted to the air is great enough to ionize it; (2) metal is ablated 
from the missile; (3) the missile is in the ionosphere where the air is 
ionized by radiation from the sun; or ( 4 ) other special devices are used 
to ionize the air, such as irradiation by high energy particles. 

The electrical conductivity of air behind a normal shock wave has 
been measured in a shock tube and is reported an page 339 of reference 3. 
The results, replotted in sketch. (z), indicate conductivities of the order 
of 1 to 5 mhos/inch for air behind strong shock waves. 

Firings of ultraspeed pellets, references 35 through 38, indicate 
that another phenomenon will enter the problem of conductivity. Photo- 
graphs and spectrographs taken with special cameras show that metal is 
eroded or ablated from the pellet and swept off the model into the wake 
where it bums. Since the metal is in a liquid spray or vapor form while 
it is in the vicinity of the pellet, the conductivity there would be 
increased by a sizable amount. It would be difficult to estimate the 
orders of magnitude to be expected. The tests show that the air can be 
made a conductor by ionization of its molecules and atoms and/or by metal 
vapor or spray eroded from the surface of the vehicle. It Is pointed out 
in reference 39 that the latter effect predominates when meteors enter the 
earth* s atmosphere. 

4 Good electrical conductors are always good thermal conductors. The 
ratio may be expressed approximately as (ref. 7)> 



where k is the coefficient of thermal conductivity, k Is Boltzmann's 
constant, and e is the charge on an electron. 
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The air In the Ionosphere Is partially Ionized "by radiation from the 
sun so that it is a conductor without addition of energy from the aircraft. 
The amount of ionization at various altitudes and the total number of 
particles are listed in tables 2 and 9 of reference 40 . A rougi estimate 
of the conductivity is found from reference 7 as 

n 2 7 jr 

ff = — — — » 3x10 ~ ie n — mhos/inch (67) 

rngUe N 


where I e , 1%, and U e are the mean free path length, mass, and velocity, 
respectively, of the electron and U 0 and N are the particle density at 
sea level and the altitude in question. The factor n is the electron 
density. The values for cr range from If )" 4 mhos/in. at 62 miles altitude 
to about 3 2 mhos/in. at 250 miles altitude. 

The large variety in fuels makes it difficult to bracket the magnitude 
of the conductivity of combustion products but information is available for 
several cases. It was reported in reference 5 that the conductivity of the 
Jet from a 50 -pound thrust acid-aniline rocket motor was measured as being 
about 10 “® mhos /in. It was added however, that the addition of a amaH 
amount of an impurity such as sodium would strikingly increase the electron 
density. 

A discussion is given in reference 6 of similar experiments and theory 
on the combustion of hydrocarbon fuels. The experimental evidence required 
that the number of electrons in the combustion products (lQ 11 to 10 14 
per cc) be several orders of magnitude greater than predicted theoretically 
( 10 s to 1 CP per cc). The maximum conductivity can be estimated by 
equation (67) to be roughly 10" 2 mhos/inch. 

Although the conductivities discussed are low compared with those 
of metals, they appear to be large enough so that appreciable magneto- 
hydrodynamic forces can be exerted. If the conductivity is too low for a 
specific case, it may be possible that alkaline salts or some other 
additive can be used to increase the number of free electrons and ions. 


APPLICATIONS 


The possibility of using a magnetic field to prevent a hot fluid from 
coming in contact with a wall which it would dissolve or destroy was dis- 
cussed in the introduction. The extent to which this can be accomplished 
is measured directly by the parameter m = (aB o 2 /pu 0() ). The larger the 
parameter m Is, the more effectively the "fluid suspension" process can 
be carried out. The influence of the magnetic field on the fluid Is 
expressed in equations (33), ( 44 ), ( 49 ), (53), (58), ... by the product 
mx. The distance x represents an "adjustment distance" required before 
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the fluid responds a given 
amount or adjusts to the mag- 
netic lines of force. If the 
magnetic parameter is very 
large^ the adjustment distance 
will "be very small and the 
fluid will follow the magnetic 
lines of force quite closely 
as indicated in sketch (p). 

The magnetic force restrains 
any motion across the lin es of force because the inertia forces are 
insignificant in comparison. In this limit the magnetic lines of force 
Eire the same as Invisible solid boundaries. 



If the parameter m is very small the adjustment distance x becomes 
very large and the influence of the magnetic field becomes negligible. 

The numbers quoted in the last section for the conductivity and field 
strength indicate that most of the applications will lie between the two 
extreme limits. 


Flow Over a High-Speed Vehicle 


Magneto-aerodynamic effects may be used to alter the external flow 
over a missile. Examples which illustrate the changes in the boundary 
layer will now be discussed. 

The Rfrin friction on an aircraft ml git be reduced by Judicious use 
of a magnetic field either to inhibit transition to turbulent flow or 
to alter the boundary-layer profile. It has been shown in references 13 
and 14 that the length of i«m-t na.r iron in a channel can be extended if a 
transverse or parallel magnetic field is applied. The results for the 
parallel field may perhaps be used as a rough estimate for the boundary 
layer flowing over a missile. Care should be exercised in applying the 
results of reference 1^ because the transverse magnetic field changes the 
velocity profile so that the pressure gradient in the channel and the 
shear stress at the wall are Increased. The net result is an adverse 
effect on the friction drag. An advantage, counteracting the increased 
friction, may be realized through a stabilization of the boundary 
layer, thus avoiding transition. 

The numerical results of table I and the equations (33) and (Ml-) 
may be used to estimate the effects of a magnetic field attached to an 
aircraft on boundary layers of uniform electrical conductivity. The 
model analyzed theoretically will probably never be duplicated in a 
physical Situation. However, circumstances may arise wherein the con- 
ductivity and magnetic field are such that they vary only slightly or 
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are almo st uniform in a layer of finite depth. This might occur, for 
example, (l) in the boundary layer of high speed, aircraft; (2) in the 
region between the surface of a blunt body and. its shock wave; ( 3 ) in 
local regions near an aircraft when material is being ablated, from the 
surface; or (4) in the ionosphere wherein air is ionized, by radiation 
freon the sun (ref. 40) and. is then a uniform conductor, provided the 
influence of the aircraft is small. The fact that the layer is of finite 
thickness causes the equations to be in error at the interface. However, 
if the conducting layer is relatively think the equations would serve as 
a good estimate in the vicinity of the surface. 

Thg second semi -infinite flat plate case which was analyzed, 
equations (49) and (54), could be thought to correspond to flight over 
the magnetic poles of the earth or through some other magnetic field. 

Hie electrically conducting fluid could come about by the same processes 
mentioned in the previous paragraph. As before, the magnetic field 
strength and electrical conductivity will probably not correspond exactly 
to the model analy zed. A close approximation is obtained when a missile 
flies through the ionosphere over one of the magnetic poles of the earth 
(provided the influence of the aircraft on the electrical conductivity of 
the air and on the magnetic field is small). Since the earth’s magnetic 
field is HTnall (b q a 4 lines/sq in. from ref. 4l), effects of this type 
are generally negligible in the lower reaches of the ionosphere. At the 
higher altitudes the magnetic parameter is large enough to cause a c han ge 
in the heat-transfer rate but at those altitudes the m e an free path is so 
large that the continuum flow analysis of this report is no longer 
applicable. 

If the air in the boundary layer is 
ionized by frictional heating, the con- 
ductivity of the air will change with, 
distance from the surface (sketch (q)) . 

Hie actual variation to be expected would 
be difficult to predict theoretically. 

The solution for a boundary layer with 
variable electrical c on d u ctivity found 
in a previous section of this paper 
(eqs. ( 58 ), (59), ( 60 ), and (64)) and 
tabulated in table ITT can be used as an 
approximate ans wer. Equation ( 58 ) shows that the skin friction can be 
reduced by a sizable amount at the expense of increased total drag (see 
eq. ( 60 )). The advantages of using a magnetic field lie in the possi- 
bility of reducing the aerodynamic heating, as expressed by equation (64). 

r Php. preceding discussion dealt with the advantages that may be 
achieved by reducing the skin friction, heat transfer, or turbulence. 
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Favorable effects may also be obtained 
by Increasing the displacement thick- 
ness of the boundary layer. Such a 
device, sketch (r), could, be used, on 
the lower surface of a wing or body to 
induce forces for lift, control pur- 
poses, or drag. The magnitude of these 
forces Is easily changed by Increasing 
or decreasing the field strength which 
In turn Increases or decreases the 
displacement thickness. 


Several additional problems arise when magneto -hydrodynamic effects 
Eire used. The system generating the magnetic field will need to be cooled 
since its parts will dissipate power. Also, the weight of the unit 
required to generate a magnetic field of sufficient strength may be too 
large to consider, or of such a size that any advantages are nullified. 

If the slow moving fluid is hot enough, to radiate energy, a sizable quan- 
tity of heat may be transferred even though- the hot fluid is not in con- 
tact with the surface. An evaluation of this effect would require a more 
complete analysis than the solutions developed in this paper. 


Whether the applications Just discussed can actually be made depends 
on the magnitude of the quantities entering the parameter m. Consider, 
for example, an aircraft moving through the upper atmosphere at 25,000 
feet per second. From sketch (l), a conductivity of not over 5 mhos/inch 
la to be expected. If p/p Q is taken as 10" 4 , and the above values are 

inserted into equation (65b), 


m = 6.43X10" 10 ^2- ss IO^Bq 2 

2.5 ■ 

In sketches (f) through (o), it is seen that the product mx should be 
of the order of 0.1 for a notable magneto-aero dynamic effect to be brought 
about. If x Is taken as about 100 inches, the magnetic field strength, 
Bq, should be around 1000 lines per Bquare inch so that the product mx 
is equal to 0.1. A magnetic field of such a strength is not unreasonable . 


Flow in Propulsion Units 


The discussion on the reduction of skin friction and heat-transfer 
rate to the walls of an aircraft by electromagnetic forces applies as 
well to the internal flow of reaction engines. It may be added that a 
uniform magnetic field across the Jet stream corresponds to the channel 
flow problems studied in references 8 to 12 and Ik. In general, imposing 
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such, a field increases the skin friction, equation (k9a), and lovers the 
thrust because of the interference of the magnetic field. A compensating 
influence is realized through a reduction in turbulence in the stream 
(see, e.g. , refs. 9 and lA) -which lowers the skin friction, heat transfer, 
and probably the noise level of jet engines. Shallow magnetic fields 
appear to be the roost promising for reducing the heat transfer if the flow 
is nearly laminar. 

The cases and illustrations considered so far have assumed that an 
electric potential E has not been applied across the fluid. To do so 
would add another force term (see eq. (2)) to equations ( 7 ) and (ll). If 
a source of power is available, a combined electric and magnetic field can 
be used to maneuver the fluid. A possible application would be to use the 
combined fields to induce an electric motor action to obt ain very high 
exhaust velocities. 


CONCLUDING- REMARKS 


The laminar boundary-layer solutions presented in this paper for the 
flow of an electrically conducting fluid over a flat plate indicate the 
changes that will be brought about by a transverse magnetic field. It 
has been found that the skin friction and heat-transfer rate are reduced 
if there is no relative motion between the plate and the magnetic field. 
The skin friction and heat-transfer rate increased in the case where 
relative motion was assumed. Irf all the cases studied the total drag was 
increased. 

The product of the magnetic parameter m and an adjustment distance 
x is the quantity which determines idle effectiveness of the magnetic 
field to chan ge the heat-transfer rate. Other studies have shown that 
a magnetic field will inhibit transition to turbulent flow and thereby 
retain a lower heat-transfer rate. 

All the cases studied, except one, assume that the conductivity and 
m ag ne tic field, did not change with distance from the surface. A more 
realistic case of a variable conductivity was considered in one example 
where a special variation was chosen. The reduction in heat-transfer 
rate and skin friction for a given value of mx was not as large as 
when the conductivity was constant. This points oxrb the importance of 
considering the way in which the conductivity and magnetic field strength 
change with distance from the surface. 


Arnes Aeronautical Laboratory 

National Advisory Consul ttee for Aeronautics 
Moffett Held, Calif., Mar. 13 , 1957 
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(a) First-order stream 
function 
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2.000 ' 

.001 

2.002 

-.002 

2.000 

.001 

2.001 

-.001 

2.000 

0 

2.000 

-.001 

2.000 

0 

2.000 

-.001 

2.000 

0 

2.000 

-.001 

2.000 

0 

2.000 

-.001 

2.000 

0 

2.000 

-.001 
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TABLE H. - MAGNETIC FIELD FIXED RELATIVE TO FLUID 


(a) First-order stream 
function 

(l) Temperature function 
ftor cold plate 

(c) Temperature function 
for Insulated plate 

m 





Ta* 

*a 

t 2 * 

D 

fa 

fa 

fa 

^ 2 /2Cp 

Uoo 2 /2Cp 

^ 2 /2Cp 

«oo a /2Cp 

0 

0 

0 

1.147 

0 

2.022 

3.228 

0 

.2 

.022 

.210 

.953 

.337 

1.370 

3.163 

-.631 

.4 

.081 

.382 

• 773 

• 555 

.821 

2.985 

-1.121 

.6 

. 17 ? 

.520 

.606 

.672 

.370 

2.723 

- 1.477 

.8 

.287 

.626 

•453 

.709 

.014 

2.403 

-1.710 

1.0 

.420 

.702 

.312 

.684 

-.253 

2.047 

-1.830 

1.2 

.566 

.752 

.186 

.613 

-.438 

1.677 

- 1.851 

1.4 

.719 

• 777 

.072 

• 513 

-.549 

1.332 

-1.786 

1.6 

.876 

! .781 

-.028 

• 398 

-596 

.967 

-1.652 

1.8 

1.031 

.767 

-.113 

.279 

-.589 

.654 

- 1.466 

2.0 

1.181 

.737 

-.184 

• 175 

-.541 

.383 

- 1.245 

2.2 

1.325 

.694 

-.241 

.064 

-.464 

.158 

-1.006 

2.4 

1.458 

.642 

-.282 

-.019 

-.369 

-.019 

-.766 

2.6 

1.581 

.582 

-.310 

-.083 

-.269 

-.149 

-.539 

2.8 

1.691 

.519 

-.323 

-.127 

-.172 

-.237 

-.337 

3-0 

1.788 

.454 

-.324 

-.153 

-.087 

-.286 

-.166 

3-2 

1.873 

• 390 

-.314 

-.163 

-.017 

-.305 

-.030 

3-4 

1.944 

.329 

-.296 

-.161 

.036 

-.301 

.069 

3-6 

2.004 

.272 

-.270 

-.149 

.072 

-.280 

.134 

3.8 

2.054 

.221 

-.240 

-.133 

.092 

-.249 

.171 

4.0 

2.093 

.176 

-.208 

-.114 

.099 

-.213 

.184 ' 


2.124 

.138 

-.176 

-.094 

.097 

-.177 

.1B1 

ESI 

2.149 

.105 

-.145 

-.075 

.089 

-.142 

.166 

139 

2.167 

.079 

-.117 

-.058 

.077 

-.111 

.144 

4.8 

2.181 

.058 

-.092 

-.044 

.064 

-.084 

.120 

5.0 

2.191 

.042 

-.071 

-.033 

.051 

-.063 

.097 

5.2 

2.198 

8 

0 

-053 

-.024 

.o 4 o 

-.045 

.076 

5.4 

2.203 

.021 

-.039 

-.017 

.030 

-.032 

.057 

5-6 

2.206 

.014 

-.028 

-.012 

.022 

-.022 

.042 

5.8 

2.209 

.009 

-.020 

-.008 

.016 

-.015 

.030 

6.0 

2.210 

.006 

-.013 

-.005 

.011 

-.010 

.021 

6.2 

2.211 

.004 

-.009 

-.003 

.007 

-.006 

.014 

6.4 

2.212 

.002 

-.006 

-.002 

.005 

-.004 

.010 

6.6 

2.212 

.001 

-.004 

-.002 

.003 

-.003 

.006 

6.8 

2.212 

.001 

-.002 

-.001 

.002 

-.002 

.004 

7.0 

2.213 

0 

-.002 

0 

.001 

-.001 

.002 

7-2 

2.213 

0 

-.001 

0 

.001 

-.001 

.002 

7-4 

2.213 

0 

-.001 

0 

0 

0 

.001 

7-6 

2.213 

0 

0 

0 

0 

0 

.001 

7.8 

2.213 

0 

0 

0 

0 

0 

0 

8.0 

2.213 

0 

0 

0 

0 

0 

0 

8.2 

2.213 

0 

0 

0 

0 

0 

0 

8.4 

2.213 

0 

0 

0 

0 

0 

0 

8.6 

2.213 

0 

0 

0 

0 

0 

0 
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TABIE HI.- MAGNETIC FIELD FIXED RELATIVE TO ELATE - 
VARIABLE CONDUCTIVITY 


(a) First-order stream 
function 



0 

-.006 

-.025 

-.055 

-.095 

-.145 

-.202 

-.264 

-.330 

-.398 

-.466 

-.532 

-.596 

-.655 

-.709 

-.757 

-.799 

-.835 

-.866 

-.891 

-.911 

-.927 

-.940 

-.949 

-.956 

-.961 

-.965 

-.968 

-.970 

-.971 

-.972 

-.972 

-.972 

-.973 

-.973 

-.973 

-.973 

-.973 

-.973 

-.973 

-.973 

-.973 

-.973 

-.973 


0 

-.062 

-.122 

-177 

-.226 

-.267 

-.300 

-.323 

-.336 

-.341 

-.337 

-.325 

-.307 

-.283 

-.256 

-.226 

-.196 

-.167 

-.139 

-.113 

-.090 

-.071 

-.054 

-.041 

-.030 

-.022 

-.015 

-.011 

-.007 

-.005 

-.003 

-.002 

-.001 

-.001 

-.001 

0 

0 

0 

0 

0 


-0.313 

-.307 

-.289 

-.262 

-.226 

-.185 

-.139 

-.092 

-.045 

0 

.004 

.077 

.106 

.128 

.143 

.150 

.150 

.144 

.134 

.121 

.106 

.090 

.075 

.060 

.047 

.036 

.027 

.020 

.014 

.010 

.007 

.005 

.003 

.002 

.001 

.001 

0 

0 

0 

0 

0 

0 

0 

0 


(l) Temperature function 

for cold plate 

1 

Te'.:. { 

WESiSmWZZSM 

0 

-0.103 

-.012 

-.021 

-.009 

.056 

.009 

.122 

.039 

.173 

.077 

.207 

.121 

.223 

.165 

.220 

.208 

.199 

.244 

.163 

.272 

.117 

.290 

.064 

.298 

.010 

.294 

-.041 

.281 

-.085 

.261 

-.118 

.235 

-.140 

.205 

-.151 

.175 

-.152 

.146 

-.145 

.118 

-.132 

.093 

-.115 

.072 

-.097 

.054 

-.079 

.040 

-.063 

.029 

-.048 

.020 

-.036 

.014 

-.027 

.009 

-.019 

.00 6 

-.013 

.004 

-.009 

.003 

-.006 

.002 

-.004 

.001 

-.003 

.001 

-.002 | 

0 

-.001 

0 

-.001 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


(c) Temperature function 
for Insulated plate 



NACA - Lui*l«y Field, Vt 
















